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We introduce the idea of effective dark matter halo catalog in f{R) gravity, which is huilt using the effective 
density field. Using a suite of high resolution N-body simulations, we find that the dynamical properties of halos, 
such as the distribution of density, velocity dispersion, specific angular momentum and spin, in the effective 
catalog of f[R) gravity closely mimic those in the ACDM model. Thus, when using effective halos, an f(R) 
model can be viewed as a ACDM model. This effective catalog therefore provides a convenient way for studying 
the baryonic physics, the galaxy halo occupation distribution and even semi-analytical galaxy formation in f{R) 
cosmologies. 


Introduction. It has become well established that the Uni¬ 
verse is currently undergoing a period of accelerated expan¬ 
sion null. The predominant explanation for this phenomenon 
is that it is driven by a non-zero cosmological constant, A, 
in the framework of General Relativity (GR). Together with 
the assumption that most of the matter in the Universe is cold 
and dark (non baryonic), this forms the current ACDM stan¬ 
dard cosmological model. There are however theoretical argu¬ 
ments as to why such an explanation should be disfavoured, 
such as the discrepancy between the value of the cosmolog¬ 
ical constant measured astronomically and that predicted by 
quantum field theory (see, e.g., Ref. Qa for a review). An 
alternative explanation is that GR might not be accurate on 
cosmological scales and so some modification to it may be 
necessary to match observations. 

A popular family of modified gravity models come under 
the umbrella of chameleon /(i?) gravity (so called because 
it replaces the Ricci scalar R in the Einstein-Hilbert action 
in GR with some function /(i?); see, e.g., Ref. ifTTl for 
recent reviews). This function introduces an effective cosmo¬ 
logical constant, which allows the Universe to expand in a way 
to match observations |[I3][T1, as well as an extra scalar de¬ 
gree of freedom, which mediates a fifth force. However, it also 
contains an efficient screening mechanism which can suppress 
this fifth force in high density environments ESlIIhl, therefore 
mimicking GR in environments such as our solar system and 
the early universe. 

Placing constraints on such modified theories of gravity can 
help us understand a lot about the nature of the cosmic accel¬ 
eration. This is one of the main tasks of upcoming cosmolog¬ 
ical surveys such as the Euclid mission IITtII . In order to make 
competitive forecasts for constraining f{R) with upcoming 
surveys we need mock galaxy catalogs which resemble the ob¬ 
servations we expect the satellite to make. These are needed 
to investigate systematic errors which could impact observa¬ 
tions. To do this, we must not only eventually be able produce 
large N-body simulations covering many Gpc^h~^ but also to 
populate these dark matter simulations with galaxies. 

In order to produce mock galaxy catalogs in modified cos¬ 
mologies, it is necessary to have an understanding of galaxy 
formation and evolution. Although galaxies are extremely 


complicated objects and many details of the physical pro¬ 
cesses still remain poorly understood even within the ACDM 
paradigm, encouraging progress has been made in recent 
years. State-of-the-art hydrodynamical simulations such as 
ILLUSTRIS m and Eagle 09], with proper modelling of 
subgrid astrophysics, are able to reproduce galaxy properties 
that are in good agreement with observations. Although hy- 
drodynamical simulations can faithfully follow the “gastro- 
physics” in a gravitational field during the hierarchical process 
of structure formation, such simulations are computationally 
expensive and high resolution simulations in f{R) gravity are 
not available currently. An alternative approach is to use semi- 
analytical galaxy formation models ||20l, which derive galaxy 
properties from dark matter simulations. Processes such as 
the cooling of gas, the formation of stars, feedback effects, and 
galaxy mergers closely relate to the properties of their host ha¬ 
los (e.g., halos mass, velocity dispersion and merger history). 
Although semi-analytical models usually contain free parame¬ 
ters, their predictions are found to be in reasonable agreement 
with observations and the models are well motivated by un¬ 
derlining physics. It is therefore of particular interest to study 
semi-analytical galaxy formation models in f{R) gravity. 

However, modifications to gravity increase the complexity 
of galaxy formation. In f{R) gravity the properties of halos 
depend not only on their mass but also on their level of screen¬ 
ing. For instance, the velocity dispersion is radically differ¬ 
ent in unscreened halos in f{R) gravity compared to halos 
of equivalent mass in ACDM. Consequently, the virial tem¬ 
perature of gas in these halos is higher than in ACDM 121]. 
It therefore follows that a halo which has assembled under 
enhanced gravitational forces may have altered astrophysics 
and it should not be assumed that the processes which govern 
galaxy formation are the same in such a halo. 

In order to overcome these difficulties, in this paper we in¬ 
troduce the idea of the effective halo catalog, which is built us¬ 
ing the effective density field in f{R) gravity. We shall show 
that the dynamical properties of halos in this catalog closely 
resemble those in ACDM dark matter halos. 

Setup. The formation of large-scale structure in f{R) grav¬ 
ity is governed by the modified Poisson equation, 


^ 2 . I&ttG , 

V , 


Email address: jianhua.he@brera.inaf.it 


3 


( 1 ) 




2 


as well as an equation for the scalar field fa, 

^^SfR= ^[SR-SirGSp] , (2) 

where (jj is the gravitational potential, Sfn = fniR) — fR{R), 
6R = R — R, and 5p = p — p. The overbar denotes the back¬ 
ground values of quantities, and V is the derivative with re¬ 
spect to physical coordinates. Combining Eq. Q and Eq. Q, 
it follows that 

= AttGSp , (3) 


where 


4>l = 4> + 


C^SfR 
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is the lensing potential. The gravitational potential (j) is felt by 
massive particles and is therefore the potential associated with 
the dynamical properties of halos and the processes of galaxy 
formation in f{R) gravity. 

We ran a suite of high-resolution N-body simulations using 
the ECOSMOG code E^ . itself based on the publicly available 
N-body code RAMSES fTij . to solve Eqs. •EH- We studied 
an f{R) model which exactly reproduces the ACDM back¬ 
ground expansion history m. Our simulations have a box 
size of Lbox = 64/i“^Mpc and contain N — 256^ particles. 
The background cosmology matches the Planck lO best-fit 
ACDM model (0° = 0.049,0° = 0.267,0° = 0.684, h = 
0.671, ris = 0.962, and erg = 0.834). Initial conditions, at a 
redshift of z = 49, were generated using the MPGRAFIC pack¬ 
age Il24l . Eourteen simulations were run in total, one realisa¬ 
tion for f{R) models with Jrq = — 10“® and Jrq = — 10“®, 
and five for fRo = —10“®'® (where /rq is the present value 
of d//di?). Eor each f{R) simulation we ran a ACDM one 
with the same initial conditions as a control. 

Effective halo catalog. We define an effective density field 
Spefi SO that the modified Poisson equation, Eq. 0, in /(i?) 
gravity can be cast into the same form as Eq. 0 

V2</> = inGSpeS , (4) 


where Spes = (| — jSUSp^^P ■ identify halos in simula¬ 
tions using the true density field Sp and the effective density 
field 6peB, respectively. The halo radius i?h is defined as the 
radius of a sphere within which the average density, is 
times the mean density, pm- The total mass inside the halos is 


Mb = 



a well defined lensing mass Ml = / 5p(x)dV and dynamical 
mass Mr = J Speff(x)dV. However, there is not a one-to- 
one correspondence between the masses in the two halo cat¬ 
alogs and there are three aspects of differences between the 
two catalogs: given the same halo radius, the lensing mass in 
the standard catalog is slightly different from that in the effec¬ 
tive catalog; the positions of centers of halos in the effective 
catalog are different from those in the standard catalog; the 
number counts of halos in the two halo catalogs are totally 
different as well. Eurther, it should be noted that Ml and Mr 
as defined above are not dependent on the shape of the halo. 
The dynamical mass Mr can be calculated accurately by our 
definition without the approximation that halos are spherical 
(e.g., in Refs.ESma). 

Scaling relations. After defining the standard and effective 
catalogs, we investigate the relationship between the mass and 
velocity dispersion, cr^, of halos in these two catalogs. The 
virial temperature of gas in virilized gaseous halos is related 
to the velocity dispersion via a power law, which also applies 
to f{R) gravity 11211 . We can thus infer the virial temperature 
of gas in dark matter halos by studying the M-a^ relationship. 
The 3D velocity dispersion of a dark matter halo in the halo- 
rest frame is defined by 


2 _ 




where Vh and u) are the halo and particle velocities, respec¬ 
tively, and Np is the number of particles inside the halo. 

Eigure shows the M-a^ relation for f{R) models with 
= -10-6, _l0-5-5,-10-5 at z = 0. In the left-hand 
panels, the mass used is the lensing mass Mr from the stan¬ 
dard catalog; in the right-hand panels, the mass is the dynam¬ 
ical mass Mr from the effective catalog. The points represent 
/(i?) halos and the color indicates their level of screening, 
with the ratio Mr/Mr illustrated in the color bar to the right. 
The black crosses represent the halos in the ACDM simula¬ 
tions, and the red and black dashed lines represent the mean 
values. In the standard catalog (left-hand panels), we can see 
that the velocity dispersions of the well-screened halos (blue) 
overlap with ACDM halos of equivalent mass. Eor unscreened 
halos (red), which are in general less massive, the M-a/, rela¬ 
tionship is different, with a yjd/'i enhancement in the velocity 
dispersion compared to the ACDM case. However, when we 
plot the velocity dispersion against the dynamical mass in the 
effective catalog (right-hand panels), the M-a/, relationship is 
the same as in ACDM for all halos. 

Next, we investigate the relationship between halo mass and 
angular momentum, which is defined by 


We modified the Amiga Halo Einder (Ahf) ll25l to identify 
dark matter halos and remove unbound particles taking into 
account the modification of gravity. Throughout this work we 
take Ah = 328 and limit our study to halos containing more 
than 400 particles. We call the catalog of halos identified us¬ 
ing the true density field 5p the standard catalog. In contrast, 
we call the catalog of halos identified using the effective den¬ 
sity field (5peff the effective catalog. The standard and effective 
catalogs are two different catalogs. In each catalog, a halo has 


Np 

J = ^ mi AT* X Aff , (5) 

i 

where Aff and Aff are the position and velocity of the z-th 
particle relative to the mean value of the halo, rrii is the true 
{inertial) mass of the z-th particle. The magnitude of the an¬ 
gular momentum is defined as J = (Jj -f -f In the 

literature, people often use the specific angular momentum, j 
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FIG. 1. The scaling relation of velocity dispersions with respect to halo mass. The points represent f{R) models. The color represents the 
ratio Mb/Ml- The black crosses represent the ACDM model. The red and black dashed lines are the averaged values. In the left panels, 
velocity dispersions in f{R) gravity do not scale as a power-law with the lensing mass. The scaling also depends on the screening. In the right 
panels, the scalings in f{R) gravity are the same as that in the ACDM model. 


= J/M, where M is the inertial mass of halos. Therefore we 
use the lensing mass Ml in the definition of j in f{R) gravity. 

The upper panels of Fig. show the scaling relation of the 
specific angular momentum j relative to the mass of halos for 
the f{R) model with fuQ = —10“^'® at z = 0. The triangles 
represent the mean values in each mass bin. In the standard 
catalog, /(i?) gravity boosts the mean value of j by a factor of 
a/4/3 for completely unscreened halos (M < ||28l, 

though this is far outweighed by the large scatters around the 
mean. In the effective catalog, f{R) gravity slightly lowers the 
mean value of j with a small relative difference Aj/jacdm ~ 
3% within the mass bin [10^^-®, 

In the literature, the spin parameter is usually defined 
asE9l. 


V2VR 


( 6 ) 


where j is the specific angular momentum and V is the circu¬ 
lar velocity at radius R — R^- The lower panels of Figure]^ 
show ln(A) as a function of the halo mass. In ACDM, A de¬ 
pends weakly on the halo mass IMI . Within the mass bin 
[10^^-®, the mean value of A is very small, A « 

0.0336 (compared to, say, a self-gravitating and rotationally- 


supported disk in which A « 0.4) with a relatively large scat¬ 
ter, (Tin A ~ 0.617. The smallness of A indicates that dark 
matter halos are mainly supported by random motions of their 
particles rather than by coherent rotation. 

In f{R) gravity, A also weakly depends on the halo mass. 
In the standard catalog, we use the definition of A as 


MLff2GMLR 


(7) 


to be consistent with the one used in the literature(e.g. 
Ref ESlO . The averaged spin parameter is significantly 
boosted by the fifth force as A « 0.0367, with a rela¬ 
tive difference AA/Aacdm ~ 9.2% within the mass bin 
[10“/1013-5]Mq. 

In the effective catalog, the circular velocity V at radius R 
is defined using the dynamical mass as = GMd/R (see 
Eq. (|4|). The spin parameter A, therefore, 

A = .- . (8) 

Mls/2GMdR 


The average spin parameter is very close to the ACDM value, 
A « 0.0326, with a small relative difference of AA/Aacdm ~ 
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FIG. 2. Upper panels: the scaling relation of halo mass with the spe¬ 
cific angular momentum j = J/Ml- The black triangles represent 
the mean values in each mass bin for the ACDM halos. The blue and 
red triangles represent the mean values in each mass bin for f{R) 
halos in the standard and effective catalogs, respectively. The dashed 
lines indicate the la scatter. Lower panels: the logarithm of the spin 
parameter ln(A) as a function of the halos mass. 

3%. Both standard and effective catalogs have approximately 
the same size of scatter, at « 0.617. 

Profiles. We now turn to the halo profiles of density, ve¬ 
locity dispersion, and specific angular momentum. Again, we 
focus our study to the f{R) model with /^o = —10“^'® at 
z = 0. We choose two different mass bins in which most dark 
matter halos are unscreened, 0.95 — 1.05 x {10^^, IO^^IMq 
( see Fig.[^. We only consider profiles at r > 10 kpc//i since 
the accuracy of the halo profiles below this radius is affected 
by the limited resolution of our simulations. 

The velocity dispersion as a function of the radius is defined 
as 

P ieAr 

in which ANp is the number of particles within the spherical 
shell Ar at a given radius r. 

The top panels of Fig. show the density profiles for ha¬ 
los in the two mass bins. The solid lines represent the mean 
values and the shaded regions represent the Icr scatter. We 
can see that the density profiles of halos in f{R) gravity and 
the ACDM model are indistinguishable in Icr range of scat¬ 
ter. On the other hand, when we look at the velocity disper¬ 
sion profiles (middle panels), we can see that in the standard 
catalog there is an enhancement in f{R) gravity compared 


FIG. 3. Top panels: the density profiles for different mass bins. 
Middle panels: the velocity dispersion profiles. Bottom panels: 
j = JjMcir < i?) within a sphere as the function of radius. The 
solid lines are the averaged values. The shaded regions represent the 
la scatter. 


to ACDM. However, we can see that the velocity dispersion 
profiles are almost identical in both the effective and ACDM 
catalogs. Finally, in the bottom panels of Fig. we show the 
specific angular momentum j(< r) as a function of distance 
from the halo center. The halos from the effective catalog and 
the ACDM model show a good agreement in the distribution 
of j. 

Halo mass function. In Fig. we show the mass function 
of the effective and standard halos compared to that of the 
ACDM model. It is clear that the effective halos show more 
significant enhancement in their abundance than that of the 
standard halos. This result indicates that the statistics of stan¬ 
dard halos may seriously underestimate the impact of f{R) 
gravity on galaxy clustering. 

Summary. In the standard catalog of f{R) gravity, halo 
properties depend not only on their masses but also on their 
level of screening. However, by introducing the effective dark 
matter halo catalog, we find that the relationships between the 
effective mass of a halo and its dynamical properties, such as 
the density profile (or equivalently potential profile), velocity 
dispersion, specific angular momentum and spin, closely re¬ 
semble those in ACDM cosmology. Thus, the effective dark 
matter halo catalog provides a convenient way for analyzing 
the baryonic physics in f{R) gravity such as the X-ray clus¬ 
ter and the Large-Scale Sunyaev-Zel’dovich effect lISTl [^ . 
Since the effective halos show more significant enhancement 
in their abundance than that of the standard halos, the mea- 
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FIG. 4. Mass function of the effective and standard halos compared 
with that of the ACDM model. ninM = is the comoving 

number density of halos per logarithmic interval in mass M. The ef¬ 
fective halos show more significant enhancement in abundance than 
that of the standard halos. 


surement of galaxy cluster abundance can provide more ro¬ 
bust constraints on f{R) gravity than using the standard halos 
(e.g. Ref. Il33l .) The effective dark matter halo catalog can 


also give us some basic insights into the galaxy formation in 
/(i?) gravity. On the scale of galaxies, the self-gravity of gas 
(which is the dominant baryonic component of galaxies) can 
be neglected, and so “gastrophysics" such as cooling and the 
accretion of gas, in an effective halo in f{R) gravity should 
follow the same relationship with mass as halos in a ACDM 
cosmology. It can therefore be expected that galaxy halo occu¬ 
pation distribution models, designed to work in a ACDM cos¬ 
mology, can be straightforwardly applied to the effective halo 
catalog in f{R) gravity. Furthermore, although we demon¬ 
strated this idea for a specific f{R) gravity model, it can be 
generalised to other modified gravity or coupled dark energy 
models, and therefore is expected to have a much wider appli¬ 
cation. 

Of course, the baryonic physics inside galaxies can also be 
changed by modified gravity, especially in regions where the 
self-gravity of baryons dominate over the dark matter. The 
processes such as the formation of stars and feedback might be 
sensitive to the modification of gravity, and as a result galaxy 
properties such as color and luminosity may differ from the 
ACDM predictions. Although further studies on these topics 
are needed, we have seen that the introduction of effective halo 
catalog can greatly simplify the analysis of physical processes 
of galaxy formation in modified gravity and therefore will be 
a useful initial step in this direction. 
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